Impurity states in multiband s-wave superconductors: analysis of iron pnictides 
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We examine the effect of a single, non-magnetic impurity in a multiband, extended s-wave super- 
conductor allowing for anisotropy of the gaps on the Fermi surfaces. We derive analytic expressions 
for the Green's functions in the continuum and analyse the conditions for the existence of sharp 
impurity-induced resonant states. Underlying band structure is more relevant for the multiband 
than for single band case, and mismatch between the bands generically makes the formation of the 
impurity states less likely in the physical regime of parameters. We confirm these conclusions by 
numerically solving the impurity problem in a tight-binding parameterization of the bands relevant 
to pnictide superconductors. 



PACS numbers: 

Since the discovery of the multiband iron-based su- 
perconductors (FeSC) 1-6 , the determination of their gap 
shape has been an active subject of study. Several theo- 
ries proposed that the superconducting order parameter 
for many pnictides and dichalcogenides is of the s-wave 
(Ai g ) type. However the presence of several Fermi sur- 
face (FS) sheets away from the center (r-point) of the 
Brillouin Zone (BZ) allows for strongly anisotropic, or 
nodal (with zeroes) gap. Experiments sensitive to low 
energy quasiparticle excitations indicate that the gap 
structure is non-universal across the families and dop- 
ing ranges, so that complementary methods have to be 
employed for a complete picture to emerge. 

These measurements include the temperature depen- 
dence of the penetration depth and spin-lattice relax- 
ation rate, both temperature and field variation of the 
specific heat and the thermal conductivity, and angle- 
resolved photocmission, among others. Scanning tunnel- 
ing spectroscopy (STS) at or near impurity sites played 
an important role in testing the anisotropic supercon- 
ducting state of the high-T c cuprates 7 , and therefore it 
is natural to ask what information it can provide for the 
iron-based materials. In conventional single-band super- 
conductors with an isotropic gap, potential scatterers do 
not change the local density of states appreciably. In 
contrast, for superconductors with sign-changing order 
parameter on the Fermi surface, non-magnetic impuri- 
ties create a quasi-bound (resonant) state, whose energy 
relative to the gap maximum and the shape in real space 
both carry information about the gap shape. 

Existence of multiple FS sheets complicates the pic- 
ture, and several recent studies arrived at different con- 
clusions. For tight-binding bands with the extended s- 
wave cos(fc x ) cos(fc y ) gap (in the "unfolded" BZ, see be- 
low), Refs. 8,9 do not find low-energy impurity states, 
and Rcf. 10 finds such states at or above the energies of 
about half of the SC gap amplitude. In two-band con- 
tinuum models of isotropic gaps of opposite signs (s^ ), 

Refs. 11,12 find that the impurity resonances can form 
deep within the gap. The Bogoliubov-de Gennes analy- 



sis of the five-band model 13 finds the low energy impu- 
rity states for some intermediate values of the scattering 
potential. The inevitable question is how general the re- 
sults obtained using a specific set of assumptions are, and 
whether different results are due to details of the models 
or salient features of pairing. 

We address this issue in the current Communication. 
We combine analytical and numerical techniques to in- 
vestigate the impurity resonance states for different gap 
shapes in the A\ g representation. One of our conclusions 
is that, for similar gap shapes, details of the band struc- 
ture, even for similar Fermi surface topologies, strongly 
affect the location of the resonance state. Also, the role 
of inter- and intra-band impurity scattering potentials 
in the formation of the impurity resonance is different 
depending on the degree of anisotropy in the gap. We 
explain the physics behind these effects by combining an- 
alytical and numerical approaches. 

The simplest model of FeSC superconductors that cap- 
tures the relevant physics has two FS sheets: one hole- 
like (h) around r-point and one electron-like (e) close to 
M and equivalent points in the BZ. It allows the anal- 
ysis of the two principal scenarios for the A\ g pairing: 
a) s± where magnetically-assisted predominantly inter- 
band pair scattering requires a sign change in the order 
parameter between FS sheets of different type 14 , and b) 
,s ++ where orbital fluctuations promote pairing with the 
same sign on the FS sheets 15-17 . Under both pictures the 
gap on the hole sheet is nearly isotropic in the x-y plane, 
while Coulomb repulsion may lead to an anisotropic, or 
even nodal, gap on the electron sheet. 

To describe scattering by a non-magnetic impurity we 
make a simplifying assumption that the scattering am- 
plitude depends only on the band index of the initial 
and final states, so that H imp = J^kk'a Ujfc' jk(T Cj>k<a , 
where Ujji = Uo if the band indices j = f (j = e, h) 
and Ujj> = U\ otherwise. The approximate indepen- 
dence of the elements of the 4 x 4-matrix U (we use U 
to denote 2x2 matrices in Nambu space) is justified by 
the small size of the FS in FeSC. This parametrization 
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means separation into small and large momentum trans- 
fer scattering, and naively one can expect U\ -C Uq be- 
cause of screening. However band structure calculations 
show that the same Fe d-orbitals contribute significantly 
to both the electron and the hole sheets of the Fermi sur- 
face, 18 and hence an impurity at or near the Fe site, will 
produce a significant inter-band scattering component. 
One candidate is the Co-dopants in the 122 series 19 . 

We compute the Green's function which is a matrix 
in both band and particle-hole space, Gjj> (k, k'; r) — 



det(A) = D(lu) = 0, where 



-(Tr 



*ik(r)®*t, k ,(0) 



where the Nambu spinor 



is *] k 



=( c jkt' 



Cj-ki), and T r is the imaginary time or- 
dering operator. In this notation the Hamiltonian for 
a pure superconductor is -ffjk=£jkT3 + AjiA, with To 
the identity matrix and % (i=l. . .3) the Pauli matri- 
ces in the Nambu space, £jk is the quasiparticle en- 
ergy in band j, and Ajt is the superconducting gap 
function on the j-th Fermi surface sheet. We ig- 
nore the weak dispersion along the c-axis. In the ab- 
sence of impurities Gojj> (k, k^—Sjj'S^'GojCk), with 
Go.j(k; iu> n )~(iuj n TQ — Hj±) , and Matsubara frequen- 
cies are ui n = 2nT{n + 1/2). 

For an extended s-wave Ak = A + 
B[cos(k x a) + cos(k y a)], in the "unfolded" zone scheme, 
with a is the lattice constant for the square Fe lat- 
tice. 14,18 ' 20,21 . When projected on the hole, Sh, and 
electron, S e , FS sheets this results in a nearly isotropic 
A/jk ~ A/, for k € Shi an d a generally anisotropic gap on 
the electron sheet(s), A e k=— A e (l + r cos 20). Here <f> is 
the angle as measured from the [100] and [010] directions 
at (±7r,0) and (0,±7r), respectively. In our notations 
So = Ag/A/i > for s± state and So < for the s ++ 
state. Below we compare the values r = (isotropic gap 
on electron Fermi surface), r = 0.7 (nodeless anisotropic 
gap), and r = 1.3 (nodal gap). 

We solve the single impurity problem using the T- 
matrix (T) method, where T and G satisfy the coupled 
equations 7 G = G + GoTG and f = U + U J2k &oT. 
For the momentum-independent U, the solution T = 

[I — UJ2k "ol U is solely a function of the band in- 
dex and the frequency. Upon analytic continuation 
iu n — > ui + i0 + the poles of G give the energies of 
elementary excitations, and hence the poles of T(uS) 
give the energies of the impurity-induced states. The 
density of states per spin in each band is Nj(r,uj) = 
— 7r _1 Im (Gjj,ii(r, r; ui + i + )) , where indices 11 refers 
to the particle component in the Nambu space. 

Denoting gj = Sk G o.j(k) = J2i9jin, with i = 
0, . . . , 3, we find that the components of the T-matrix 
satisfy T ee = Uqt 3 + U T 3 g e T ee + Uyr 3 g h T he and T he = 
U1T3 + Ui? 3 g e T ee + UoT 3 g h T he ; the equations for T hh 
and T e h are obtained by switching indices. The solu- 
tion is f ee = A' 1 [Uo? 3 - (U^ - U 2 )T 3 g h r 3 ] where A = 
(1 - U ? 3 g h ) (1 - U ? 3 g e ) - U 2 T 3 g h ? 3 g e . Consequently, 
the energies of the bound state are determined from 
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5J 2 =E,(-1)\9 



and 

Note that Eq. (1) de- 



we introduced 
cr-inNjUo)- 2 (1 - U g ]3 ) 2 
pends only on U 2 , and hence we take U\ > without 
loss of generality 

As in single-band superconductors, the gap shape af- 
fects the angular averages of the anomalous Green's func- 
tions, <7ji, and therefore influences the energy of the 
impurity state both in a single- and multi-band cases. 
Two other aspects reflect the multiband nature of the 
system. First, since summation over momenta in the 
T-matrix equations yields a prefactor of the density of 
states (DOS) at the appropriate FS sheet, the ratio of 
the DOS n = N e /Nh controls the number of states avail- 
able for inter- vs. intra-band scattering. In the mod- 
els with dominant intcrband pairing the same parameter 
controls the ratio of the gap amplitudes, So- In analogy 
with Rcfs. 22,23 we find that at T = such models obey 
the constraint uSq(1 + r 2 /2) = 1. Variation of n very 
significantly affects the properties of the impurity bound 
state. 

Second, the particle-hole band asymmetry, 



9j3 = 1^2 

k 



UJ: 



(2) 



appears in combination with the the impurity potential in 
the equation for the T-matrix. This term is largely deter- 
mined by the normal state band structure, and therefore 
can be approximated by its value with A^k = 0. In ana- 
lytical approaches it is often assumed gj 3 = 0, although 
even in single-band superconductors the resonance state 
is sensitive to the band structure 24 ' 25 . In FcSC the situ- 
ation is even more complex since the chemical potential 
is close to the top/bottom of the hole/electron bands re- 
spectively. As a result: a) g e3 and gh 3 have opposite 
signs and hence the sign of Uq matters; b) we expect 
9ez/N e ~ —ghz/Nh ~ 0(1), which changes the bound 
state properties relative to the particle-hole symmetric 
case. 

The location of the bound state for the isotropic gap, 
r = 0, is given by 



a 2 (A e + A^) 2 -2(a 2 -l)A e A^ 
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2(a 2 - 1) 
g(A e + A h ) 
2(a 2 - 1) 



(3) 



^V« 2 (A e - A /i ) 2 + 4A e A ?1 



where all the information about the DOS and scattering 
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FIG. 1: (Color online) Energy of the impurity bound state for 
isotropic gaps (r = 0). Here ljm = min(w+, see Eq. (3), 
and minimal gap Am = min(A e , Ah). For a given value of the 
DOS ratio N e /Nu, and the corresponding fixed <5o, variation 
of the impurity potentials Uo and Ui traces a path along the 
surface. For a given value of a, the bound state energy is 
always lowest at 5o=l. 
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Fig. 1 shows that a) the bound state is much deeper in 

the gap for the s_| than for the s++ state; b) the lowest 

energy of the bound state is reached for the "balanced" 
band case, N e = Nh and = A e . and any deviation 
from this regime leads to the impurity state edging closer 
to the continuum. Fig. 2 shows several traces over this 
surface for fixed n, where the resonance state becomes 
closer to the gap edge and exists in a narrower and nar- 
rower range of parameters as the densities of states in 
the two bands "detune" . Hence not only the topology of 
the Fermi surfaces, but also the curvature of the bands at 
the Fermi level and the relative bandwidth matter. In our 
view it is for that reason that the simplest tight-binding 
parameterization 26 does not give a clearly defined bound 
state. For that model we estimate N e (0) sw 4.9^(0), 
g e 3 « -3.77riV e (0), and g h3 « 0AirNh(0), which suggests 
that the impurity state is very close to the gap edge, or 
does not form at all, in agreement with the conclusions 
of Rcfs. 8,9. In the extreme case of the proximity to van 
Hove singularity for one of the bands, the impurity states 
form only for U\ 3> Uo 27 . In contrast, the studies utiliz- 
ing more realistic band structure fits with smaller ratio 
of the DOS often find the deeper and sharper impurity 
resonance 12 . 

For r — the lowest energy of the impurity resonance 
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FIG. 2: (Color online) Energy of the impurity bound state 
for r = for different DOS ratios, n, and the particle-hole 
anisotropics, gjz, measured in units of 7riVj. We denoted 
N F = (N e + N h )/2. Note that the deep and well-defined 
bound state only appears for the strong interband scattering 
Vi. 
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For particle-hole symmetric case (gr e 3 = gnz = 0) this 
clearly indicates U\ > Uo, which is unlikely. For strongly 
particle-hole asymmetric bands, and for strong scattering 
(c e ,/i < 1), the impurity bound state may be formed for 
more physical values of U\ < Uq, but, as is seen from 
Fig. 1, the state itself is at finite energy. 

As the value of r, and the anisotropy of the gap on the 
electron sheet increases this trend remains for as long as 
the gap does not develop nodes. Once the gap on the 
electron sheet develops zeroes, a more familiar aspect of 
the single-band nodal superconductors comes into play: 
intraband scattering produces a bound state, while the 
interband scattering broadens it. This is illustrated in 
Fig. 3. Even in this case the small size of the Fermi sur- 
faces, i.e. moderately large and opposite in sign values 
of gjz lead to a much shallower resonance state that for 
a particle-hole symmetric case. Note also a very signifi- 
cant anisotropy with respect to the sign of the intraband 
scattering potential. These analytical results show that 
in multiband systems with opposite nature of the carri- 
ers in the two bands it is generically difficult to realize a 
well-defined impurity resonance states. 

To verify that these conclusions remain valid for 
pnictides we numerically solved the T-matrix equation 
and determined the local density of states on the 
impurity site for a tight-binding fit to the Fermi surface 
from Ref. 28. In the folded (2-Fe) BZ the energies of 
the hole and the electron bands are given by £, ai k = 
— t ai (cos (k x ) + cos (k y )) — t' a . cos (k x ) cos (k y ) — fx ai , and 

C/3,k = -tfr (COS (kx) + COS (ky))-fp t COS (^f) COS J - 

fip i respectively, with the hoppings and the band 
shifts (in eV) (t ai ,t' ai ,fi ai ) = (-0.3,-0.24,0.6) 
and (t a2 ,t' a2 ,fi a2 ) = (-0.2,-0.24,0.4) for the hole, 
and (tfai/p^Ufa) = (-1.14,-0.74,-1.70) and 
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FIG. 3: (Color online) Energy of the impurity bound state 
for r = 1.3 for different DOS ratios, n, and the particle-hole 
anisotropics, gjs, in the same notations as in previous figure. 



(tfat'p^fjLfo) = (-1.14,0.64,-1.70) for the elec- 
tron bands. In this model n ~ 0.4, close to the n = 0.5 
considered above, and different from n w 0.2 used 
in Refs. 8,9, even though the Fermi surface topology 
is similar. Hence the difference between the results 
stems from this DOS imbalance. Our Green's functions 
g e3 -7riV e (0), and g h3 « 1.2^^(0). We again 
consider the isotropic gap A Qi k=A/ l on the hole FS 
sheets, and we define the gap on the electron sheets 

to be 2 k = — A e ( 1 ± f cos (^-) cos f^fjjj where 
the upper and lower signs arise from the folding of the 
Fermi surfaces into the smaller 2-Fc Brillouin Zone. In 
general the gap anisotropy on the electron FS sheets 
is not a simple cosine, but is close to it in shape. We 
can extract the effective anisotropy from the ratio 
|A e , m i„|/|A eimoa ;| = |r-l|/|r+l|, and make comparison 
with our analytical results. The electron Fermi surfaces 
in this description are very close to each other, so 
that the values of r differ by about 6% between them, 
and we quote the average number. The calculations 
were performed on a 2000x2000 /c-space grid with the 
intrinsic broadening 7=Ah/40=0.0015 eV. Fig. 4 shows 
the LDOS at the impurity site for moderate values of 
the scattering potentials. As before, the particle-hole 
band asymmetry causes the impurity-state energy to 
depend strongly on the sign of Uq. For these moderate 
values of Uq and U\ the low-energy states do not form 
when the FS sheets are fully-gapped (r = 0) except for 
the unphysical case U\ > Uq. The impurity resonance 
forms more easily when there are nodes on the electron 
sheets; this happens regardless of the sign of So since 
the formation of the state is dominated by intraband 
processes. For anisotropic nodeless gap (r = 0.7) there 
exist broad features associated with the transfer of 
spectral weight from the coherence peak to energies 
above the threshold e m = A e \r — 1|. These states mix 
with the continuum and are not sharp. At the same 
time we find that they can feature relatively prominently 
at the impurity site simply because the peak at e m is 
much smaller than that at u = Ah in the pure system. 
A sharp bound state close to mid-gap only appears for 
extremely high values of Uq ~ XJ\ ~ 100, also supporting 




FIG. 4: (Color online) Evolution of on-site LDOS (eW 1 ) 
for four-band model with moderate scattering, [/o=±1.5, for 
(a,b) isotropic, (c,d) nodeless anisotropic, and (e,f) nodal gaps 
for So > 0. Panels (g,h) are for nodal gap with So < 0. 
The thin black line is the DOS of clean system, and arrows 
mark the DOS feature at e m = (0.38, 0.41) for the anisotropic 
cases with r = (0.7, 1.4), respectively. Low-energy impurity 
states form below e m in the nodal system even at small values 
of Ui, but impurity states do not form near lj=0 in either 
of the fully-gapped systems. LDOS for Ui t m is shown when 
Ui.m 9^ \Uo\ but is unlikely when Ui t m > \Uo\- The insets 
show close ups of the low-intensity positive-bias peaks. Note 
the different vertical scales. 



the qualitative analysis above. 

Our main conclusions therefore are that within a model 
of relatively momentum-independent inter- and intra- 
band scattering in multiband systems a) the mismatch 
of the densities of states and the gap values on different 
Fermi surface sheets affects very significantly the energy 
of the impurity bound states, pinning them relatively 
close to the continuum states over most of the param- 
eter range; b) this conclusion remains qualitatively valid 
even for nodal gaps on one of the Fermi surface sheets 
for realistic particle-hole anisotropies in the electron and 
hole bands; c) broad features due to impurities may ex- 
ist for anisotropic gaps in the s± case, but are far less 
likely for the s++ pairing. We explain the differences in 
the results between different groups as stemming in part 
from the different underlying band structure. Of course, 
in a fully microscopic theories starting from the orbital 
representation the resulting effective scattering may be 
anisotropic within each Fermi surface sheet: strong be- 
tween parts with similar orbital content and weak be- 
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tween regions stemming from different Fe orbitals. It 
would therefore be very instructive to check such effec- 
tive potential, for example, in the model of Rcf.13. It 
will also be useful to check whether within these models 
the scattering potential varies strongly on the scale of the 
bandwidth, since such an effect assists the formation of 
the bound state even for an isotropic gap. Our results 
clearly show that the study of impurity states alone is not 
sufficient to draw reliable conclusions about the shape of 



the superconducting gap in multiband systems. 
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